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Abstract 

The expansion in multipoles Ji,i = 2, . . . of the gravitational po- 
tential of a rotating body affects the orbital motion of a test particle 
orbiting it with long-term perturbations both at a classical and at 
a relativistic level. In this preliminary sensitivity analysis, we show 
that, for the first time, the J2C~^ effects could be measured by the 
ongoing Juno mission in the gravitational field of Jupiter during its 
yearlong science phase (10 November 2016-5 October 2017) thanks to 
its high eccentricity (e = 0.947) and to the huge oblateness of Jupiter 
(J2 = 1.47 X 10~^). The semi-major axis a and the perijove uj of Juno 
are expected to be shifted by Aa < 700 — 900 m and Aw < 50 — 60 mil- 
liarcseconds (mas), respectively, over 1 — 2 year. A numerical analysis 
shows also that the expected J2C~^ range-rate signal for Juno should be 
as large as w 280 microns per second (/im s^^) during a typical 6 h pass 
at its closest approach. Independent analyses previously performed by 
other researchers about the measurability of the Lense-Thirring effect 
showed that the radio science apparatus of Juno should reach an ac- 
curacy in Doppler range-rate measurements of « 1 — 5 /im s^^ over 
such passes. The range-rate signature of the classical even zonal per- 
turbations is different from the IPN one. Thus, further investigations, 
based on covariance analyses of simulated Doppler data and dedicated 
parameters estimation, are worth of further consideration. It turns out 
that the J2C~^ effects cannot be responsible of the flyby anomaly in the 
gravitational field of the Earth. A dedicated spacecraft in a 6678 km x 
57103 km polar orbit would experience a geocentric J2C~^ range-rate 
shift of sa 0.4 mm s^^. 

PACS: 04.80.-y; 04.80.Cc; 95.55.Pe 
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1 Introduction 



The general theory of relativity is one of the fundamental pillars of our 
knowledge of the natural world, being, to date, the best theory of the grav- 
itational interaction at our disposal. Thus, it is important to put it on the 
test in as much ways as possible to increase our confidence in it, especially in 
view of the extrapolations of its validity to extreme scenarios. In these cases, 
empirical checks are difficult and/or the interpretation of existing observa- 
tions heavily rely upon more or less speculative assumptions concerning the 
history of the systems considered and the physics governing them in such 
regimes. 

In this respect, the weak-field and slow- motion approximation of general 
relativity played so far an important role. If on the one hand, the magnitude 
of its predicted effects is generally modest with respect to those expected in 
the strong-field regime, on the other hand, they occur in the gravitational 
fields of stars and planets such as our Sun, the Earth and some of the 
major bodies of our Solar System which can be probed with great confidence 
with man-made objects tracked with increasing accuracy. Moreover, the 
key physical parameters of such natural laboratories are, in general, well 
known from a variety of independent space-based missions, thus reducing 
the untested (and often untestable) speculations down to a minimum. 

So far, the weak- field tests of general relativity probed the behaviour of 
the electromagnetic waves propagating in the vicinity of massive objects, 
and the orbital dynamics of test particles and gyroscopes moving around 
them [l]. As far as the latter class of effects are concerned, at the first 
Post-Newtonian (IPN) level they are 

• The precession of the pericenter of a test particle in the field of a 
spherically symmetric, static source. It corresponds to the celebrated 
Einstein precession of the perihelion of Mercury [2] of 42.98 arcseconds 
per century [3]. 

• The geodestic or de Sitter-Fokker precession of a gyroscope orbiting a 
spherically symmetric, static body [IHS]. It was measured in the field 
of the Sun by tracking the motion of the Earth-Moon system with 
the Lunar Laser Ranging (LLR) technique [7H9], and in the field of 
the Earth with the man-made gyroscopes carried onboard the Gravity 
Probe B (GP-B) spacecraft HO]. 

• The gravitomagnetic effects occurring in the vicinity of a stationary 
source: the Lense-Thirring orbital precessions of a test particle and 
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the Pugh-Schiff precession of a gyroscope [T2l[T3] orbiting a rotating 
body. A measurement of the Pugh-Schiff effect was recently made 
with the GP-B mission |10] . For a recent overview of the attempts to 
measure the Lense-Thirring effect in the Solar System, see, e.g., [H] 
and references therein. 

One of the major sources of systematic bias in the attempts to detect such 
relativistic orbital effects is represented by the competing classical secular 
precessions induced by the expansion in multipoles of the Newtonian grav- 
itational potential of the central body accounting for its departures from 
spherical symmetry. It turns out that the asphericity of the gravitational 
field of the primary has an impact on the motion of a test particle also at the 
IPN level |15H20j . inducing non- vanishing long-term orbital perturbations 
on all the Keplerian orbital elements |21H25j . So far, such effects did not re- 
ceive much attention since it was always believed that they are too small to 
be detectable in any foreseeable future. For example, Brumberg [21] wrote 
that "these effects are not observable as yet and [. . . ] can hardly have a 
great chance of being detected in the near future" . 

The scope of this paper is to show that recent developments in space 
sciences may soon overturn such opinions. After reviewing the IPN orbital 
effects due to the oblateness of the central body in Section [21 the possibility 
of detecting them in some astronomical systems is considered in Section [3l 
In Section 13.11 we will show that there are good perspective to accurately 
measure the J2C~^ orbital effects in the gravitational field of Jupiter with 
the Juno spacecraft currently en route to it. Also the Earth scenario is 
examined in Section [3. 2[ Section [J] summarizes our findings. 

2 The IPN orbital effects of order J2C~^ 

The direct Post-Newtonian perturbations due to the quadrupole moment 
J2 of the central source of mass M and equatorial radius R were explicitly 
computed in alternative metric theories of gravity by Soffel et al0 [21] and by 
Brumberg |24j by means of the standard Gauss equations for the variation 

^The rates of change of eq. (|l])-eq. ([5]) can be obtained by taking the time derivatives of 
the long-periodic (LP) and secular (S) shifts Aa, Ae, . . . of the orbital elements calculated 
by Soffel et al. [2T] and by using the general relativistic values of the parameters entering 
them. 

■^For a better comparison with the results by Brumberg, the following replacements are 
to be made in his formulas at pag. 114-115 of [H]: Q -> -J2GMR^ = -J2nla'^ R ,m ^ 
GM/(?,ot — >■ 0, — j> a;. Here c is the speed of light in vacuum, G is the Newtonian 
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of the Keplerian orbital elements which arqj the semi-major axis a, the 
eccentricity e, the inclination / of the orbital plane to the reference {x, y} 
plane, the longitude of the ascending node /2 and the argument of pericenter 
u. In general relativity, they are, to first order, 

da _ 9 J2n3i?2e2 _^ ^^^2 j 2oj 

di~ 8c2 (1 - e^f ' 



de 21J2nlR^e{2 + e^)sm^ Ism2uj 
di ~ 16c2 (1 - e^f ' 

dl _ 3J2nlR^e'^ sm2I sin2uj 



dt 8c2 (1 - e^f 

dQ _ 3J2nlR^ cos / (6 - cos 2uj) 
'dt ~ 4c2 (1 _ e2)3 ' 

^ ^ 3J2nlR^ 132 _ 3e2 _ 2 (7 + 2e') cos 2u+ 

dt 32c2(l-e2)^^ ^ ^ 



(3) 



(4) 



+ cos 21 [48 - 9e^ + (l4 - 4e^) cos 2uj] } , (5) 

We performed our own calculations within general relativity by using the 
Gauss equations and the true anomaly / as fast variable of integration; by 
using a Keplerian ellipse as unperturbed reference trajectory, we were able 
to obtain eq. ([TJ-eq. ([5]). 

The indirect, second order mixed perturbations induced by the New- 
tonian quadrupole field and the Schwarzschild acceleration were computed 
in general relativity by Heimberger et al. j23j by using the Lie transform 
perturbation theory [26] with the Delaunay variables / = Ai,L = n\^a?',g = 

constant of gravitation, and Wb = V GMa~^ is the Keplerian mean motion of the test 
particle. 

^In the following we will not consider the mean anomaly M since it is rarely considered 
in empirical studies. It is so because the Keplerian mean motion should be subtracted from 
its rate; actually, the knowledge of ?ib is limited by the uncertainty in the primary's GM 
and, in the case of artificial spacecrafts, by any non-gravitational perturbations plaguing 
the semi-major axis o. 
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u),G = nba^Vl — e^, h = fi, H = ni^a^y/l — cos /. In terms of the usual 
Keplerian orbital elements, they are [23] 



da 



(6) 



sin^ / sin 2uj 
dt ~ 4c2 (1 _ e2)2 ' 



(7) 



dl 



(8) 



3J2nlR'^ cos / (^7 + 296^ + ISVl - - Se^ cos 2a; 
^ ~ 4c2 (1 - e2)3 



(9) 



^ {876^ +3(2 + e2)cos2a;+ 



dt 16c2 (1 - e2)3 



+ cos2/ 1456^ - 3 (2 + Se^) cos 2uj + 108\/l - e2 + 107 + 



+ 60^1 - e2 + 45} . 



(10) 



For another computation based on the Lie transform perturbation theory 
and the Delaunay elements, see [25] . 

Both eq. ([I])-eq. ([5]) and eq. ©-eq. pOj) hold in a frame whose reference 
z axis is aligned with the unit vector k of the spin axis of the central body. 



3 Perspectives for a detection 

In this Section we look for some astronomical scenarios which may allow 
for a measurement of the J2C~2 orbital effects. Contrary to what one might 
think at first glance, such a perspective is not unrealistic. 
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3.1 Jupiter and Juno 

Junc0 [27] is a spacecraft en route to Jupiter where its arrival is scheduled 
for July 2016. Its science phase has a nominal duration of almost 1 yr (10 
November 2016-5 October 2017) [25]. During it, Juno will move along a 
highly ehiptical (e = 0.947), polar (I = 90 defl orbit to accurately map, 
among the other things, the gravitational field of Jupiter through the 
perturbations on the Juno's wide (a = 20.03i?, = 11.07 d) trajectory 
due the departures from spherical symmetry of the Jovian gravitational 
potential. 

In view of the high diurnal rotation rate of Jupiter and of the peculiar 
orbital geometry of Juno, it was suggested [30] that it could also carry 
out a measurement of the Lense-Thirring effect at a percent level yielding 
a determination of the angular momentum of Jupiter; further independent 
analyses confirmed such a possibility |28p31| . From a practical point of view, 
the measurement would be conducted with the Juno's radio science systerqj 
(X-Band, Ka-Band) providing an accurate determination of the Doppler 
shift [28l[3l] during manjjil of the scheduled 31 orbits in the science phase, 
mainly at the perijove passages lasting about 6 hr. Indeed, as pointed out 
in [31j , large longitude-keeping maneuvers would compromise the dynamical 
coherence of the orbit, making, thus, practically infeasible an analysis based 
on steady time series of the Keplerian orbital elements. The expected overall 
range-rate accuracy should be of the order of 1 — 5 ^m s~^ over time scales of 
10^ s [28l[30l[3l] for the Ka-Band apparatus, while the spacecraft's position 
and velocity should be known with an uncertainty of about 10 m and 1 mm 
s~^, respectively, in the three spatial directions [28]. 

The high Jupiter's oblateness (J2 = 1.47 x 10~^), combined with the 
relevant eccentricity of the Juno's orbit and with its accurate radio science 
Ka-Band system, may offer a unique opportunity to detect, for the first time, 
the J2C~^ orbital effects. Here we will offer just a preliminary sensitivity 
analysis; in view of its promising outcome, it is hoped that it will prompt 

*See http://www.nasa.gov/mission_pages/juno/main/index.html on the WEB. 

^In a coordinate system whose reference z axis is aligned with the Jupiter's spin axis. 

®While the Ka-Band system is for gravity science, the X-Band apparatus is used for 
the spacecraft orbit determination and navigation [28]. As far as the Doppler range-rate 
measurements are concerned, the X-Band system is less accurate than the Ka-Band one 
by a factor of about 10 - 100 [28] . 

^It depends on the pointing of the Juno's high-gain antenna [28]; if it does not point to 
the Earth, as in the microwave radiometry passes when it will point to Jupiter, Doppler 
tracking is impossible. According to Helled et al. [28], about 22 useful perijove passes 
should be available during the scheduled science phase. 
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further, more detailed investigations by independent teams in such a way it 
occurred for pHj lST]. 

As a first outlook, in Figure [1] we depict the shifts Aa and Aw of the 
semi-major axis and of the perijove of Juno as functions of the initial value 
ujQ of the perijove itself and of the duration T of the science phase. We 
obtained Aa and Alo by integrating eq. ([1]) and the sum of eq. ^ and 
eq. (flOl) from t = to t = T. We posed to = ujq + Cjt, where oj is the sum 
of both the classical J2 and of the Schwarzschild+Lense-Thirring perijove 
precessions. Although the nominal duration of the science phase of Juno is 
about 1 year, we considered also an extension of it, as it is not infrequent 
in successful space missions. The shift of the semi-major axis can reach 
values as large as Aa ~ 600 — 900 m, while the perijove shift can amount to 
I Awl ~ 50—60 milliarcseconds (mas); from such a point of view, an extension 
of the science phase up to ~ 1.5 — 2 yr would be certainly desirable. We did 
not show the eccentricity shift Ae since it turns out to be as little as w 5 
mas; the node and the inclination are not affected at J2C~^ level because 
of the polar orbital configuration. It is important to note that, contrary to 
the perijove, the semi-major axis does experience neither long-periodic nor 
secular competing perturbations of gravitational origin; thus, there would 
be no further gravitational aliasing bias on a. As it can be inferred from eq. 
([5]) and eq. pO|) . the perijove has both a secular and a long-periodic J2C~^ 
component which would help in separating it from the biasing classical and 
relativistic (Schwarzschild+Lense-Thirring) secular precessions; according 
to |30j . the mismodelling in the Newtonian perijove precessions due to the 
even zonals of the Jupiter's gravitational potential would be as little as ~ 1 
mas yr~^. 

In order to make a closer connection with the actual direct observable, 
in the right panel of Figure [2] we plot a numerically integrated Earth- Juno 
range-rate signal A/j due to the J2C~^ acceleration over a 6 hr pass centered 
about the first perijove passage after the beginning of the science phase. 
The mathematical model of the J2C~^ acceleration used in our numerical 
integration of the equations of motion of Juno has no restrictions on the 
spatial orientation of k. This is a useful feature since it allows to use different 
coordinate systems. It is likely that the real data will be actually processed 
in the ICRF frame; in it, the Jupiter's spin axis is kx = —0.0146021, /c^ = 
—0.430337, kz = 0.90255, as it can be retrieved from the celestial coordinates 
of its north pole of rotation at epoch J2000.0 [32] . The resulting J2C~^ range- 
rate signal amounts to about 280 fim s~^, thus supporting the feasibility of 
our proposed measurement. From the left panel of Figure O we also note 
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Juno: semimajor axis shift of order 0(J2C ^) 
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Juno: pericenter sliift of order 0(J2C ^) 
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Figure 1: Left column, upper panel: plot of the semi-major axis sliift Aa, 
in m, of Juno as a function of the initial value ujq of the perijove and of the 
duration T, in yr, of the science phase of the mission. Left column, middle 
panel: plot of the semi-major axis shift Aa, in m, of Juno as a function of 
the initial value ojq for T = 1.6 yr. Left column, lower panel: plot of the 
semi-major axis shift Aa, in m, of Juno as a function of the duration T of the 
science phase of the mission for ujq = 60 deg. Right column, upper panel: 
plot of the perijove shift Au of Juno, in mas, as a function of the initial 
value cjQ of the perijove itself and of the duration T, in yr, of the science 
phase of the mission. Right column, middle panel: plot of the perijove shift 
Auj of Juno, in mas, as a function of the initial value ujq for T = 1.6 yr. 
Right column, lower panel: plot of the perijove shift Aw of Juno, in mas, as 
a function of the duration T of the science phase of the mission for ojq = 60 
deg. ^ 




Figure 2: From the left to the right: Earth-Juno range, in m, and range- 
rate, in fim s~^, due to the J2C~^ acceleration. Both were calculated by 
numerically integrating the equations of motion of the Earth, Jupiter and 
Juno with and without the J2C~^ term over 1 yr in a coordinate system 
with the ICRF equator as reference {x, y} plane. In it, the Jupiter' spin 
axis is = -0.0146021, fcj^ = -0.430337, = 0.90255, as it can be in- 
ferred from the ICRF equatorial coordinates at epoch J2000.0 of its north 
pole of rotation [32]. Both the integrations shared the same initial con- 
ditions corresponding to the beginning of the science phase (10 Novem- 
ber 2016) [28]; for Jupiter and the Earth they were retrieved from the 
WEB interface HORIZONS by NASA/JPL. For Juno we adopted xq = 
-0.0444405 R,yo = -1.3097 R,zo = 2.74684 R,xo = 19.4951 km s-^,yo = 
— 11.9024 kms~^,io = 23.5485 km s~^ corresponding to an initial polar 
orbital configuration with Wq = 5.7 deg with respect to the Jupiter's equa- 
tor [28] . The time interval of the plots, covering 6 hr, is centered about the 
first perijove passage after the beginning of the science phase. 
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that the J2C~^ range shift Ap of Juno is as large as ~ 20 m over the same 
interval. 

We stress once again that this is just a sensitivity analysis aimed to 
explore the possibility of a detection of the effect considered. It should be 
accompanied by a dedicated full covariance study with the simulated data of 
the real Doppler range-rate measurements at the perijove passages, and by 
the estimation of dedicated solve-for parameters explicitly accounting for the 
J2C~^ acceleration in the dynamical models. Nonetheless, the encouraging 
outcome of the investigations in |28p31j concerning the Lense-Thirring effect 
makes us confident that also for the J2C~^ effect a percent determination 
with Juno may be possible. 

3.2 The Earth's scenario 

Moving to the Earth's gravitational field, it may be worth of mentioning 
that the J2C~^ effects cannot explain the flyby anomaly [33]. We performed 
a numerical analysis for the NEAR probe since it was most evident during 
the flyby of such a spacecraft occurred on 23 January 1998, when the X- 
Band Doppler apparatus detected an unexplained range-rate shift as large 
as A/5 = 13.46 lb 0.01 mm [33] • We numerically integrated the equations 
of motion of NEAR with and without the J2C~^ acceleration starting from 
the initial conditions of the 1998 flyby. The resulting J2C~^ rang-rate shift 
was as little as ~ 10~^ mm s~^. 

A dedicated mission would have more chances, at least in principle. As 
an example, we considered a spacecraft endowed with a Doppler tracking 
apparatus moving along a highly elliptical (e = 0.79053) polar orbit char- 
acterized by an altitude at the perigee of /imin = 300 km and a maximum 
geocentric distance as large as rmax = 57103.2 km. By numerically integrat- 
ing its equations of motion in a geocentric equatorial {kx = ky = 0, = 1) 
frame, it turned out that the magnitude of its J2C~^ range-rate shift, re- 
ferred to the geocenter, increases at every perigee pass in such a way that 
at, say, the 25th perigee passage it amounts to —0.4 mm s~^. 

4 Summary and conclusions 

We re-examined the IPN effects on the motion of a test particle orbiting an 
oblate central body in order to check if recent advancements in space science 
and technology make them potentially measurable. 

The answer is cautiously positive thanks to the ongoing Juno mission to 
Jupiter. Indeed, its peculiar orbital configuration, the accuracy of its radio- 
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tracking apparatus and the huge Jovian oblateness make the perspective of 
measuring the J2C~^ effects on the Juno's orbit feasible during its scheduled 
yearlong science phase in 2016-2017. Indeed, the numerically integrated 
J2C~^ range-rate shifts at the perijove passes are expected to be as large as 
« 280 /im s~^; the onboard Ka-Band Doppler system should be accurate at 
a 1 — 5 fim level during such passes. The signatures of other competing 
effects such as the classical shifts due to the even zonal harmonics of the 
Newtonian component of the aspherical gravitational potential of Jupiter 
are different, thus likely allowing for an adequate separation of the signal we 
are interested in. Most of the positive results of independent investigations 
previously made by other researchers in view of a possible measurement 
of the smaller Lense-Thirring effect with Juno may be valid for the J2C~^ 
effects as well, thus making us reasonably confident about a successful de- 
termination of them. Nonetheless, a dedicated covariance analysis implying 
simulations of Doppler measurements and parameter estimation is required 
to further confirm the promising results of our preliminary sensitivity study. 

We also looked at the Earth's gravitational field. Our numerical analysis 
for the NEAR spacecraft discards the possibility that the flyby anomaly can 
be due to the J2C~^ range-rate shift. We also investigated a hypothetical 
scenario based on the use of a dedicated spacecraft orbiting along a highly 
eccentric polar trajectory by finding that, for a perigee height of 300 km and 
an apogee geocentric distance of about 57000 km, the J2C~^ range-rate shift 
would be of the order of 0.4 mm s~^. Such a figure should be confronted 
with the typical Doppler accuracy of about 0.01 mm reached in various 
past flybys by several spacecrafts. 
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